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We explore the dynamics of the r-modes in accreting neutron stars in two ways. First, we explore 
how dissipation in the magneto-viscous boundary layer (MVBL) at the crust-core interface governs 
the damping of r-mode perturbations in the fluid interior. Two models are considered: one assuming 
an ordinary-fluid interior, the other taking the core to consist of superfluid neutrons, type II super- 
conducting protons, and normal electrons. We show, within our approximations, that no solution 
to the magnetohydrodynamic equations exists in the superfluid model when both the neutron and 
proton vortices are pinned. However, if just one species of vortex is pinned, we can find solutions. 
When the neutron vortices are pinned and the proton vortices are unpinned there is much more 
dissipation than in the ordinary-fluid model, unless the pinning is weak. When the proton vortices 
are pinned and the neutron vortices are unpinned the dissipation is comparable or slightly less than 
that for the ordinary-fluid model, even when the pinning is strong. We also find in the superfluid 
model that relatively weak radial magnetic fields ~ f 9 G (f 8 K/T) 2 greatly affect the MVBL, 
though the effects of mutual friction tend to counteract the magnetic effects. Second, we evolve 
our two models in time, accounting for accretion, and explore how the magnetic field strength, the 
r-mode saturation amplitude, and the accretion rate affect the cyclic evolution of these stars. If the 
r-modes control the spin cycles of accreting neutron stars we find that magnetic fields can affect the 
clustering of the spin frequencies of low mass x-ray binaries (LMXBs) and the fraction of these that 
are currently emitting gravitational waves. 



I. INTRODUCTION 

The r-modes are oscillation modes that occur in ro- 
tating fluids due to the Coriolis effect. Great interest in 
these modes was generated after Andersson and Fried- 
man and Morsink [0 showed that gravitational-radiation 
backreaction tends to drive these modes unstable at all 
angular velocities. However, internal dissipation most 
likely completely suppresses this instability in all stars 
except neutron stars (|,||,^,Q,||J|J^] . The r-mode insta- 
bility in neutron stars is further complicated by the fact 
that the problem splits into two cases: superfluid neutron 
stars, with interior temperatures below approximately 
10 9 K, and within which the interior of a neutron star is 
expected to contain regions of superfluid neutrons mixed 
with lower concentrations of superconducting protons, 
normal electrons and other exotic particles, and ordinary- 
fluid neutron stars with interior tempertures above 10 9 K. 
Furthermore, in a rotating neutron star the superfluid 
neutrons will form a dense array of quantized vortices, 
and if an interior magnetic field exists the superconduct- 
ing protons (for type II superconductivity) will form a 
dense array of flux tubes (also called vortices in this pa- 
per). (See ||||l0 n], p^ , [l3t for discussions of superflu- 
idity in neutron stars.) A third case, not considered in 
this paper, is that in which the nucleons dissolve into a 
soup of up, down, and strange quarks, which results in a 
strange star, not a neutron star. 

The regimes within which unstable r-modes can ex- 



ist in neutron stars have been narrowed down as the- 
oretical understanding has improved. First, Bildsten 
and Ushomirsky |l4] showed that when a solid crust is 
present, the shear dissipation in the viscous boundary 
layer (VBL) that forms at the crust-core interface greatly 
suppresses the r-mode instability. Neutron stars are ex- 
pected to form a solid crust for p < 1.5 x 10 14 g/cm 3 
and for tempertures below an approximate melting tem- 
perature of T = 10 10 K |l5|,16[ . This work was ex- 
tended by Andersson et al. |H] , Rieutord [Q , Levin and 
Ushomirsky fljl , Lindblom, Owen, and Ushomirsky [pfj[ , 
and Mendell | ]21| . Bildsten and Ushomirsky |14j also pre- 
dicted that magnetic fields would be important in the 
VBL, for a magnetic field strength B and temperature 
T, when B > 10 11 G (10 s K/T). Mendell §| confirmed 
this and showed that magnetic fields further suppress the 
instability in ordinary-fluid neutron stars. Second, Jones 
p2| , p3[ and Lindblom and Owen |24|] have shown that 
hyperon bulk viscosity further suppresses the instability 
above a temperature of 10 9 K, while Haensel, Levenfish, 
and Yakovlev Q have shown that superfluidity of the 
baryons tends to suppress hyperon bulk viscosity below 
10 9 K. In this paper we ignore hyperon bulk viscosity 
altogether. Finally, Wu, Matzner, and Arras [^6) and 
Arras et al. J27J have shown that even if the r-modes are 
driven unstable, the saturation amplitude is likely to be 
very small. 

However, the r-modes still remain important for the 
same reasons that sparked the initial interest in them. 
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First, it is still possible that t hey play a role in the 
spin-down of hot young pulsars [P8| . Second, they could 
be reponsible for the clustering of spin frequencies in- 
ferred from th e observations of low mass x-ray bina- 
ries (LMXBs) 
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|^J^.y,^,|l7[|3jj35|]. (However, note 
that some complications exist with the interpretation of 
the observations; see |56| and references therein.) Fi- 
nally, even given the recent results, it is still possible 
that gravitational-waves from unstable small amplitude 
r-modes could be detected by enhanced or advanced 
narrow-banded gravitational-wave detectors. Theoreti- 
cal understanding of the r-modes instability is far from 
complete, and thus more work is needed to understand 
these issues. 

The purpose of this paper is to explore how the 
magneto-viscous boundary layer (MVBL) that forms at 
the crust-core interface in the presense of a radial mag- 
netic field affects the dynamics of the r-modes in accret- 
ing neutron stars. Thus, this paper focuses on neutron 
stars with a solid crust (T < 10 10 K), and we assume that 
viscous dissipation in the MVBL is the dominant form of 
dissipation in these stars. Two models are considered: 
one assuming an ordinary-fluid interior; the other, which 
we refer to as the superfluid model, taking the core to 
consist of superfluid neutrons, type II superconducting 
protons, and normal electrons. First, we explore how 
dissipation in the MVBL at the crust-core interface gov- 
erns the damping of r-modc perturbations in the fluid 
interior. This extends the work of Mendell to the 
superfluid case. (Also, a minor coding error in Mendell 
pl[ caused the MVBL damping times to come out about 
53% percent too small and the critical angular velocities 
to come out about 11% too large. We give the corrected 
results in this paper.) Second, we evolve the r-mode am- 
plitude, the angular velocity, and the temperature of the 
two models in time, accounting for accretion. We explore 
how the magnetic field strength, the r-mode saturation 
amplitude, and the accretion rate affect the cyclic evo- 
lution of these stars. This extends the work of Levin 
p3| , Anderson, et al. [Jill, and Wagoner, Hennawi, and 
Liu (see also Heyl fl35|). We use the equations in 
Wagoner, Hennawi, and Liu ]34j] to evolve the models. 

We show, within our approximations, that no solution 
to the magnetohydrodynamic (MHD) equations exists in 
the superfluid model when both the neutron and proton 
vortices are pinned. However, if just one species of vor- 
tex is pinned, we can find solutions. When the neutron 
vortices are pinned and the proton vortices are unpinned 
there is much more dissipation than in the ordinary-fluid 
model, unless the pinning is weak. When the proton vor- 
tices are pinned and the neutron vortices are unpinned 
the dissipation is comparable or slightly less than that 
for the ordinary-fluid model, even when the pinning is 
strong. We also find for the superfluid model that rel- 
atively weak radial magnetic fields ~ 10 9 G (10 s K/T) 2 
greatly affect the MVBL. We find that magnetic fields 
tend to make the critical angular velocity for the onset of 
the r-mode instability temperature independent and in- 



crease the dissipation rate, as Mendell |^T) found in the 
ordinary-fluid case, though the effects of mutual friction 
tend to counteract the magnetic effects. Even when the 
magnetic field is zero, the correct scaling of the bound- 
ary layer thickness with the proton mass density is given 
here for the superfluid model for the first time. 

When we evolve our two models in time we find that 
the critical angular velocity decreases with temperature 
sufficiently for all reasonable magnetic fields to produce 
the thermal run-away found by Levin |33[ |. (This hap- 
pens even though magnetic fields tend to flatten the crit- 
ical angular velocity vs. temperature curves. Wagoner, 
Hennawi, and Liu j34j have shown that if the critical an- 
gular velocity is (mostly) temperature independent then 
the r-mode amplitude and temperature oscillate with a 
period of hundreds to thousand of years, while the spin 
frequency of the star stays roughly constant.) Even with 
thermal run-way, Levin J33] | and Anderson, et al. p"^ ] 
showed that the r-mode instability can still produce clus- 
tering of the spin frequencies in LMXBs. If true, we 
find that magnetic fields can have important effects on 
this clustering. For radial fields of B > 10 11 G in the 
ordinary-fluid model or for B ~ 10 9 G in the superfluid 
model the spin cycle of an LMXB becomes thinner, which 
would cause LMXBs with these fields to cluster into a 
narrower range of spin frequencies. Cumming, Zweibel, 
and Bildsten |}7j have shown that buried fields in the 
crusts of accreting neutron stars must be less than 10 11 G, 
while typical external fields in LMXBs are 10 8 " 9 G. Fur- 
thermore, the interior magnetic field may be expelled 
during neutron star spindown or affected by accretion 
p8| , [39 40, 4l]]. However, we know that neutron stars typ- 
ically start out with large fields ~ 10 12 G. Since mag- 
netic field evolution is uncertain, at some point it may 
be natural for the interior field to spend time at the val- 
ues given here. If the magnetic field narrows the spin 
cycle of an LMXB, this increases the fraction of time it 
spends spinning down and radiating significant amounts 
of gravitational radiation. This, in turn, has an effect on 
the number of LMXBs that are currently radiating. Low- 
ering the saturation amplitude or increasing the accretion 
rate produces the same effects, making it difficult to in- 
fer the interior field from current observations. However, 
the detection of gravitational waves from the r-modes of 
an accreting neutron star would determine the saturation 
amplitude (if the distance to source is known). In this 
case it might be possible to place limits on the interior 
magnetic fields of LMXBs with known accretion rates 
based on the observed allowed range of spin frequencies. 

The next section reviews the Newtonian MHD equa- 
tions for superfluid neutron stars. These reduce to the 
ordinary-fluid case when appropriate limits are taken. 
Sec. IH finds an approximate solution to the MHD equa- 



tions valid in the MVBL. Sec. IV presents the results of 
calculations for the MVBL damping times and the crit- 
ical angular velocity, and discusses limits on the mag- 
netic field and the MVBL length-scales. The models 
are evolved and the spin cycles of LMXBs are studied 
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in Sec. Conclusions are discussed and suggestions for 
future work are made in Sec. [v| . 



II. SUPERFLUID MHD EQUATIONS 

The Newtonian MHD equations for a mixture of su- 
pcrfiuid neutrons, type II superconducting protons, and 
normal electrons have been derived by Mendell . The 
quantities in the equations represent smooth averages 
over volumes containing many vortices, based on the ap- 
proach of Bckarevich and Khalatnikov |^] and as ex- 
tended to neutron stars by Mendell and Lindblom p4| , 
Mendell @|6), and Lindblom and Mendell It is 

convenient to define the average and relative velocities, 
v and w respectively, in terms of the mass densities and 
velocities of the neutron and protons: 



pv = p n v n + p p Vp, 
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(2.1) 



(2.2) 



where B n is the mutual friction coefficient. In Eqs. (2.3) 
and (2A) SU and 5(3 are related to the perturbed pres- 
sure, gravitational potential, and chemical potentials in 
the star but play no further role in t his pa per (see Lind- 
blom and Mendell j47|). In Eqs. (0) and (|J) 7 is a di- 
mensionless factor defined in Lindblom and Mendell ^] 
that arises due to the "entrainment" effect that occurs 
due to strong interactions between neutrons and protons 
p9| , [T0 12 1 . When entrainment occurs the mass currents 
of the neutrons and protons are given in terms of the 
superfluid velocities by, 



Mn — PnnVn PnpVp, 



A'lp — Pnp^n ~\~ PppVp- 



(2.7) 



(2.8) 



The coefficients, p nn , p np , p PP form what is called the 
mass density matrix. These are related to the ordinary 
mass densities and the entrainment factor 7 by, 



Pnn Pn Pnp , 



(2.9) 



For a star rotating uniformly in equilibrium with angular 
velocity f2 (ft ~ |f2|), the equations in Mendell (4^] for 
small Eulerian perturbations (prefixed with a S) can be 
written in the corotating frame as 



d t Sv+2fl x Sv 




— (V x SX P ) x 



m p c 



-D 



— V • (2r] e Sa e ), 
P 



dtSw+2-fn x 5w = -V5/3 + — ( — X B 

Pp\ c 



(2.3) 



1 



(V x SX P ) x 



-B 



Pp 

-— V • {2r) e 5a e ) - -^SF n 
P P Pv 



d t SB = V x (dv e X B). 



(2.4) 



(2.5) 



For simplicity, terms of order p e / ' p p have been ignored, 
and only the largest vortex force (dependent on X p ) has 
been retained. Dissipative effects due to the shear of the 
electron fluid, a e , and the mutual friction force caused 
by electron scattering off the neutron vortices, F n , are 
included. The mutual friction force is given by p6[ 



SF n = 2njB„ 



Sw 



n(n ■ Sw) 



(2.6) 



Ppp — Pp Pnp j 



and 



Pnp 



PnPp(l -7) 



(2.10) 



(2.11) 



Next, note that Easson j5(| calculates that the ratio 
of the magnetic diffusion time-scale to the viscous dif- 
fusion time-scale in a typical neutron star is roughly 
10 14 (10 8 K/T) 4 . Thus, the electrical conductivity of the 
electrons in the core and the crust is approximated as 
infinite; magnetic diffusion plays no role for the temper- 
atures and time -sca les of interest in this paper, and is 
ignored in Eq. (|J). (See Mendell Q for further diss- 
cusion of this issue.) Finally, e is the absolute value of 
the charge of the electron, m p is the mass of the proton, 
and c is the speed of light. 

The MHD limit is valid for studies of oscillations with 
phase velocities much less than the speed of light, fre- 
quencies much less than the plasma and cyclotron fre- 
quencies, and large conductivities. Under these circum- 
stances the above equations, along with the mass con- 
servation laws and equations of state, completely deter- 
mine the dynamics of system. All the other vector fields 
of interest are determined in terms of Sv, Sw, and SB. 
Specifically, for phase velocities much less than the speed 
of light, the displacement current can be ignored in Am- 
peres law, and the current density is given by 



5 J = —V x SB. 

47T 



(2.12) 



For an infinitely conducting crust, there will also be a 
surface current density at the crust-core interface, given 
by 
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SI 



c 

4tt 



SB x f 



(2.13) 



where R c is the radius of the core. For frequencies much 
less than the plasma and cyclotron frequency, the per- 
turbed charge density is neglibible, and the perturbed 
electrical current density is so small that the electron ve- 
locity is approximate ly SM v /p v fl5|,[42|]. Using Eqs. <\2.l\] 



(f2~2|) and Eqs. (p~8|) - (P ■ 1 1|) this can be written as 



Sv e = Sv + —rySw. 
P 



(2.14) 



Finally, when the above approximations hold and the 
conductivity is high, electrons (being the least massive 
charge carrier) respond to make the Lorentz force on 
them negligible, and the electric field is given by 



5E = - S -^xB-*x5B. 

c c 



(2.15) 



Turning to the vortex force in Eqs. (2.3) and ( |2.4| ), this 
force is due to the underlying array the proton vortices 
(much smaller forces due to neutron vortices are ignored 
in this paper). It is given in terms of the vector vortex 
"chemical potential" A a (basically the energy per unit 
length needed to increase the number of proton vortices 
by one in direction a). Mendell Jl5|,fl2| shows that this 
force is given by 



(V x SX p ) x — B 



SJ 



2m Ti 



V x S(s p -) 



x B 



-B 



(2.16) 



The first term on the right side of this equation is not an 
approximation. It occurs in the exact form of the proton 
vortex force, and it always cancels the Lorentz force (oc 



5 J X B) in Eqs. fl2.3| ) and fl2.4j). For this reason, Alfven 
waves do not occur in an exotic type II superconducting 
proton, normal electron plasma. Instead, these waves are 
replaced by cyclotron- vortex waves (see Mendell |42| and 
the next section) . In the second term on the right side of 
Eq. (2.16), s p is the energy per unit length of a proton 
vortex. Small corrections to this term have been ignored 
(See Mendell p^]). Finally, throughout this paper note 
that B=\B\. 



III. APPROXIMATE MAGNETO- VISCOUS 
BOUNDARY LAYER SOLUTIONS 

Approximate solutions to the superfluid MHD equa- 
tions are found in this section. The ordinary-fluid limit 
of the solutions is also found. The equilibrium magnetic 
field is assumed to be arbitrary (for now) except that it 
is static in the corotating frame, and it is restricted such 
that no equilibrium electrical currents exist. This implies 



the equilibrium structure of the star is unchanged by the 
presence of the magnetic field. To facilitate the manipu- 
lation of tensor quantites, a rotating coordinate basis will 
be used, and indices will be raised and lowered using the 
flat-space metric tensor in spherical coordinates. In this 
basis, the equilibrium velocity is v a — Q,(f) a . Following 
the notation of previous studies, note that Latin indices 
are space indices, except where it is understood that n, p, 
e, o, and c refer to neutrons, protons, electrons, ordinary, 
and the crust respectively. 

Let Sv a , Sw a , and SB a be the standard r-mode so- 
lution valid in the bulk of the core where viscous and 
magnetic forces are small compared to the Coriolis force. 
(As explained at the end of the last section, note that 
the magnetic forces are due entirely to the proton vortex 
force in the superfluid model.) By definition, the words 
"standard r-mode solution" mean the r-mode solution 
that would exist in the fluid core, if no solid crust were 
present. For the purposes of this paper, these solutions 
are also taken to be the lowest order form of the stan- 
dard solution, when the solution is expanded in powers 
of the angular velocity. The standard r-mode solution 
is already known from previous studies of the r-modes. 
Note that SB a is approximately zero, and is taken to be 
exactly zero in the standard r-mode solution. Further- 
more, Lindblom and Mendell E8[ show that Sw a — at 
lowest order, and thus the lowest order superfluid r-mode 
solution is identical to the lowest order ordinary-fluid so- 
lution described in previous r-mode papers. (In general, 
see reference [Q for the standard r-mode solution in 
the notation used in this paper.) However, because of 
the solid crust, boundary conditions must be applied to 
the tangential components of the velocities at the crust- 
core interface. An ordinary viscous fluid cannot slip at a 
perfectly rigid boundary, and superfluid vortices cannot 
move if perfectly pinned at the boundary. However, elec- 
trons can slip at a conducting boundary. When fluids or 
vortices cannot slip at a boundary this causes a boundary 
layer to form. In the boundary layer the magnitudes of 
the viscous, magnetic, and Coriolis forces become com- 
parible. Mutual friction forces can be important too. 
Thus, all these forces have effects on the structure of the 
boundary layer. When both viscous and magnetic forces 
exist we refer to the boundary layer as a MVBL. 

Let 6v a , Sw a , and SB a be the corrections that must 
be added to the standard r-mode solution to enforce the 
tangential boundary conditions at the crust-core inter- 
face. The problem in this section reduces to finding solu- 
tions for the corrective quantities, Sv a , Sw a , and SB". 
Since the total fields are 5v a + Sv a , Sw a + Sw a , and 
SB a + SB a , and the equations are linear, the correc- 
tive quantities obey the same equations as the standard 
quantities, Eqs. (2.3)-(2.5). However, approximations 
that hold true for the corrective quantities in the MVBL 
are made when solving for these quantities that differ 
from the approximation made when finding the stan- 
dard quantities. (This is done such that total solution 
is approximately valid inside and outside the MVBL.) 
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The main assumption, following the work of Lindblom, 
Owen, and Ushomirsky |2p| and references therein, is 
that the boundary conditions force the corrective quan- 
tities to rapidly change in the radial direction. Thus, it 
is assumed that terms involving radial derivatives of the 
corrective quantities dominate the angular derivatives of 
these quanties, and dominate any derivative of all other 
quantities. (For certain angles, certain caveats must be 
added to this statement. These angles are discussed in 



Sees. IV.) By construction, the ratio of radial derivative 
terms kept to the terms dropped will be the ratio of the 
core radius to the boundary-layer length-scales. Thus, 
the main assumption is valid as long as the boundary 
layer length-scales are much less that the core radius. 

Thus, it is possible to find approximate equations for 
the corrective quantities, 5v a , 8w a , and 6B a by taking the 
following steps. First, let all perturbed quantities in the 
corotating frame have time dependence exp(iK57t), where 
k is a constant that gives the mode frequency in terms 
of ft. Second, note that since in the standard solution 
6v r , Sw r , and SB r vanish everywhere, mass is conserved 
at the boundary and the divergence of the magnetic field 
is zero if the radial components, 5v r , 5w r , and SB r also 
vanish. Finally, of the spacial derivative terms, only the 
highest order radial derivitives acting on the corrective 
quantities are kept, all other spacial derivative terms are 
dropped. Taking thes e ste ps, the resulting equations are 
as follows. First, Eq. (2.5) becomes, 



SB a = 



-B r d r 6v a e . 



(3.1) 



Substituting this, Eq. Q2.16D , and Eq. fl2J) into Eq. flgj ) 
and Eq. (2.4) gives 



p(iKfl5v 9 - 2r2cos6>sin6»<5^) / p p = Td 2 8v 6 e 



p^nnS^ + 29.cot65v B )/p p = Td 2 r 5vf, 



(3.2) 



(3.3) 



where the quantum of flux is $o = hc/2e, and rj e is the 
electron viscosi ty (a nd recall that e p was previously de- 
fined after Eq. §T|). 

Equating the right sides of Eqs. (3.2) and (3.4), and 



the right sides of Eqs. (3.3) and (3.5), the components 
Sw and 8w^ are given algebraically by 

inttppSw — 2 r yflp p cos9sm9Sw^ + 2"fQpB n Sw e 

= piinflSv 6 - 2ftcos6>sin6)<K^), (3.8) 



and 



iK^ppSw* + 2-fQ,p p coW5w e + 2'ynpB n Sw' p = 

P^kCIS^ + 2flcot6Sv e ). (3.9) 



Combinin g th ese with Eq. (2.14) shows that Eqs. ( |3.2| ) 
and Eqs. ( |3.3| ) are a 4th order system for 8v e and 8vf. 

In deriving the above equations a factor of cos 2 # 
has been ignored in the mutual friction force term in 
Eq. Q. Also, 



'B*\ SB a B b SB b 



(3.10) 



which appears in the proton vortex force, has been ap- 
proximated as 8B a /B. In the final equations this corre- 
sponds to ignoring terms of the order (B r /B) 2 (B e / B 2 ) 2 , 
{B r /B) 2 {B e B*/B 2 ) : and (B r / B) 2 (B^ / B) 2 . These 
rather crude approximations make the analysis much eas- 
ier to understand. The equations and solutions without 
these approximation are given in the Appendix. (In the 
final analysis, all magnetic terms, included the ignored 
ones, vanish if B r vanishes; B r controls the magnetic ef- 
fects on the MBVL. Thus keeping only the largest terms 
that depends on B r is not too bad.) 

The corrective boundary layer solution is then found 
by allowing all perturbative quantities to vary as 
exp[ifc(i? c — r)] (recall R c is the radius of the core). It is 
easy to show that solutions exist for 



innSw 9 - 27ftcos6»sin6»<5M)* 

+ (2 1 ripB n /p p )Sw e = Td 2 T 8v b e 



(3.4) 



k ± =K ± 



9. 



\ kS! V B ) 



PP 



(3.11) 



inflSw' 1, + 2 1 VLcoi65w e 

+ ^npBJp^Sw* = Td 2 r 5vf, 



where 



(3.5) R = ( K ± 2cos6Q[« ± 27COS6I - {2i P1 B n /pp)] 

± V [(Pn7 + P P )/p]«±27Cos0-2i 7 6„ ' 1 ' ' 



^CV f Br 



kQ \ B 



Pp 



(3.6) 



In this equation, VJ? V is the square of the cyclotron- vortex 
wave speed, defined by B2| 



T/ 2 _ £ P B 



$0Pp' 



(3.7) 



and 



Sw b 



5v 6 = ±isin06v' p , 
Sw e = ±ism9dw^ . 

p(k ± 2cos0) 
P p (k ± 27COS0) — 2ip r fB. l 



-Sv u 



(3.13) 
(3.14) 

(3.15) 
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Choosing solutions where Im(fc±) > 0, so that the solu- 
tion decays exponentially as r — > 0, the general solution 
of the equations is of the form 



5v° = [C + e lk+{R ^ r) +C-e ik -( R "- r) ]e l 



(3.16) 



We choose to define C± as the constants that give 8v®. 
Once these are determined by the tangential boundary 
conditions at the crust-core interface, all other compo- 
nents and velocities are determined in terms of these con- 
stants. 

It is easy to find the analogous solution for a mixture 
of ordinary-fluid neutrons, protons, and electrons. In this 
case, all the fluids act as a single fluid, and the e quat ions 
for this case can be obtained from (3.2) and (|3.3|) by 
taking the following limits: 





>o, 


pp — 


P, 


Ve -> 


V, 






Ve~ 


v a , 


V c 2 v 





2 

A' 



(3.17) 



where V£ is the square of the Alfven wave speed, defined 
as 



V. 



4irp 



(3.18) 



This equation should be compared with Eq. (3.7) for 
Vqv Choosing solutions where Im(k±) > 0, so that the 
solution decays exponentially as r — > 0, the general solu- 
tion for the ordinary-fluid equations is of the form 



5v b 



where 



[ C+e ik° + (R c -r) +c _ 



gife° (R c — r)lgiKf2t 



(3.19) 



k° ± 



n 



\ kO V B I 



P 



and 



\J k ± 2cos#. 



(3.21) 



These equations are exactly those obtained by Mendell 
II- 

We must now consider what the correct tangential 
boundary conditions are at the crust-core interface. In 
the ordinary fluid case it is the no-slip boundary condi- 
tion for a viscous fluid 



[Sv a + Sv a = (1 - S)Sv a } r=Rc . 



(3.22) 



Following Levin and Ushomirsky JT9[ |, but introducing 
our own notation, we introduce the "slip factor" S into 
Eq. ( p2| ). This is a simplified way to account for the 
motion of the crust. When 5 = 1 the crust is perfectly 



rigid, and no slipping of the fluids at the boundary is 
allowed. A value of S = would correspond to a fluid 
crust, in which case the corrective solutions vanish. How- 
ever, a realistic crust is not a fluid and not perfectly rigid. 
Levin and Ushomirsky ]nj found for a toy model that 
0.05 < S < 1. 

Note that in the superfluid case the simple no-slip 
boundary condition does not apply. The reasons are the 
neutron and protons are superfluid and thus can slip, and 
the electrons are conducting and thus can slip along the 
conducting crust too. Instead, the tangential boundary 
conditions at the crust-core interface are determined by 
the pinning of the neutron and proton vortices. However, 
if both are pinned then the problem is over-determined 
and no solution exist. The reason for this is that in our 
approximation scheme the only independent velocity is 
that of the electrons; the other velocities are determined 
algebraically in terms of the components of 5v e . Thus, 
if both the neutron and protron vortices are pinned, the 
approximation scheme used here breaks down and either 
the MVBL becomes much more complicated or the r- 
modes are drastically changed. 

However, if just one species of vortex is pinned, we can 
find solutions. (More realistically, such solutions proba- 
bly hold if the pinning of one vortex species completely 
dominates over the other species.) Using the equations 
in the Appendix of Mendell Eq] we find, for a particu- 
lar species, that the smooth-averaged perturbed vortex 
core velocity is proportional to the perturbed mass cur- 
rent projected perpendicular to the equilibrium direction 
of the vortex array. Since the radial components of the 
velocities are zero, this amounts to applying the tangen- 
tial bounda ry co nditions to the mass currents themselves. 



Using Eqs. fl2.7| ) and (|2.8|), pinning just one vortex species 
corresponds to the pinned-vortex boundary condition, 



g v * + e^5v» + ei^Sv" = (1 - S)5v a 

Pn Pn 



(3.20) for pinned neutron vortices, and 



Pp 



Pp 



(1 - S)5v a 



r=R c 



r=R c 



(3.23) 



(3.24) 



for pinned proton vortices. We introduce the factor S 
again to account for slipping of vortices or motion of the 
crust. A value of S — 1 corresponds to strong (perfect) 
pinning while a value of S = corresponds to no pinning. 
Small values of S correspond to weak pinning. Several 
studies imply that vortex pinning in neutron stars, in 
general, may be moderate or weak, |5l],^2|,[53| though 
certain configurations may result in increased pinning 
strength Q. (Also see Ruderman, Zhu, and Chen 
for a disscussion of pinning and vortex interactions in the 
core.) 

Applying the no-slip or pinned-vortex boundary con- 
ditions give 



C± = --SA±{6v e ± isineSv 4 '), 



(3.25) 



G 



where 



Using the definition 



1. 



(3.26) 



in the ordinary fluid case or when the proton vortices are 
pinned, and 



A, 



[{jpn + P p )/p]k ± 27 cos 9 - 2ijB n 

(Pnp/Pn)K-2ijl3n 



(3.27) 



when the neutron vortices are pinned. In this latter case, 
this factor represents the fact that since equal but op- 
posite forces couple the neutrons to the charged fluids 
(via the entrainment effect and mutual friction) the ac- 
celerations and motions of these fluids are inversely pro- 
portional to their mass densities. Thus, the factor A± is 
quite large when the neutron vortices are pinned. In con- 
trast, since the electrons move with the proton mass cur- 
rents, A± is 1 when the proton vortices are pinned. Sim- 
ilarly, since all the fluids move together in the ordinary- 
fluid case A± is 1 in this case as well. Furthermore, note 



that Sv and Sv^ in Eq. (3.25) are the same for the su- 
perfluid and ordinary-fluid models because, as already 
pointed out, the lowest-order standard r-mode solutions 
are identical for these two cases. These components of 
the velocity in the corotating frame (with their time de- 
pendence cancelled) are 



Sv° = -iAr ri 



sin 



sinflfo* = Ar m - 1 f&a m - l 0coBBe im *. 



(3.28) 



(3.29) 



As shown in previous papers, the viscous damping rate 
is given by integrating the shear over the core: 



(3.30) 



where E is the energy of the mode as defined in e.g., 
Lindblom and Mendell pS| ], but limited to the core. The 
largest contribution to the integral comes from the radial 
derivatives of the corrective boundary layer velocities, so 
that in the superfluid case 



ah - \Rl{\dr5v 6 e \ 2 + |d r sin^| 2 ), (3.31) 



ignoring terms smaller than these by a factor of the 
boundary layer length-scales over the core radius. Thus, 
it can be shown that the solutions presented in this sec- 
tion give a VBL damping time of 



2tt 2 m+3 (m + l)! 
■q e S 2 l m(2ra+ 1)!! 



R, 



2m +2 



dr, (3.32) 



where 



1 = 



2TT fTT 



I A, 



2 d+(l - cos6 



|A_| 2 |/c_| 2 d_(l + cos(9) 2 sm^ m - l 6d9d(f). (3.33) 



k - 2lT + i 
A± d± ' 



(3.34) 



the boundary layer thicknesses in Eq. (3.33) and bound- 
ary layer wavelengths are defined as 



2tt 



Re(fc±)' 



Im(fc-| 



(3.35) 



(3.36) 



The results given here reduce to those given in Mendell 
pl| in the ordinary-fluid limit. 

Before we study in detail how MVBL damping affects 
the r-modes, we note the following useful analytic expres- 
sions for t v when the mutual friction coefficient B n is 0. 
For simplicity we choose B = B r = constant. This is un- 
realistic (and unphysical unless as many field lines enter 
the crust core interface as leave it) but Mendell |2l] has 
shown that this simple model gives the same qualitative 
results as a more complicated dipole field. In this case, 
the f actor s next to K± in Eq. (3.11) and next to K± in 
Eq. ( 3.20| ) are independent of 8, and we can factor them 
out of the integral X, giving us 



1 



s 2 



CoFRelJqoFB^+iQ/T 2 



+L»OFlm I y qoFB 

for the ordinary-fluid model and 
1 



• in/i? 



(3.37) 



C SF Re{ x /qsFB + in/Ti 



D SF lm[ Jqs F B + mm 



(3.38) 



for the superfluid model. In these equations qoF = 
1/(47tk?7o), and q S F = e p /($ K?7 e ,o), where 770 and 77^0 
are the temperature independent part of the viscosities 
such that ri = ?7o(10 10 K/T) 2 and i le = ?7 e ,o(10 8 K/T) 2 . 
The rest of the constants, C... and D..., reduce to inte- 
grals over the angular variables that can be done numeri- 
cally. Furthermore, when mutual friction is added to the 
superfluid model we find the following holds to about 4 
parts in 10 4 : 



A SF Tj 

s 2 n 



Re JqspB + m/Tt 



(3.39) 



We find the constant A$f (which is positive and depends 
on B n ) b y fitting this formula to the results given by 
Eq. (pj - 



7 



IV. BOUNDARY LAYER DAMPING TIMES AND 
THE CRITICAL ANGULAR VELOCITY 



(4.3) 



One goal this study is to deterine how magnetic fields 
change the cricital angular velocity needed for the on- 
set of the r-modes instability via their effect on the 
MVBL damping rate in accreting neutrons stars, such 
as those in LMXBs. Bildsten and Ushomirsky pre- 
dicted that magnetic fields would be important when 
B > 10 11 G (10 8 K/T) and this was confirmed by Mendell 
plf for the ordinarly-fluid model. However, Cumming, 
Zweibel, and Bildsten |57| have shown that buried fields 
in the crusts of LMXBs are less than 10 11 G, while typi- 
cal LMXBs have external fields that are about 10 8_9 G. 
A minor coding error in Mendell |2l[] caused the MVBL 
damping times to come out about 53% percent too small 
and the critical angular velocities to come out about 11% 
too large, without changing the qualitative results. In 
this section we present corrected results for the ordinary- 
fluid model and extend those results to the superfluid 
model. 

As shown in previous studies, gravitational radiation 
emitted by the r-modes always tends to drive these modes 
unstable. The onset of the instability occurs when the 
dissipation rate equals the gravitational-radiation growth 
rate, tgr. Thus, for our models we calculate the critical 
angular velocity for the onset of the instability by solving 



T GRi 



(4.1) 



where r„ is given by Eq. ( pi^ ) [see also Eqs. (^37|)-(^39|)1 
and tgr for the r-modes is given by 



tgr = tgr 



2m+2 



(4.2) 



where il a = y^rGp, G is the Newtonian gravitational 
constant, and d is the average density of the star (see 
references |]|,|||). 

Attention is restricted to the case used in previous 
studies: the m = 2 r-mode for a 1.4M Q n = 1 polytrope, 
as described in, e.g., Lindblom, Mendell, and Owen Q, 
and Lindblom and Mendell . Adoption of this model 
allows easy comparison with previous studies. (Note also 
that the m — 2 case is the most susceptible to the r- 
mode instability.) For this case, kq = 2/3, the stellar 
radius is 12.53 km, and il = 8413 s . Note that the 
maximum angular velocity, where mass shedding occurs 
at the equator, is roughly 2f2 Q /3 (for any equation of 
state). The density at the crust-core boundary is given 
approximately by 1.5 x 10 14 g/cm 3 (see fl5| , |l6| ). Us- 
ing this density, the characteristic gravitational-radiation 
growth time is tqr = 4.25 s |2(| , the core radius is 
R c = 11.01km, the proton density at the boundary is 
p p = 6.6 x 10 12 g/cm 3 , and the superfluid entrainment 
factor is 7 = 1.90498 (for an entrainment parameter of 
0.04, as defined in Lindblom and Mendell Q|). In gen- 
eral, the dimensionless mutual friction coefficient is given 

by mm, 



Thus, at the crust-core boundary the mutual friction co- 
efficient isi3 n = 9.5xl0~ 5 . Furthermore, the energy per 
unit length of a proton vortex is given in Mendell pfij ] 
and in a more convenient form in Mendell p2[ . However, 
even the latter form depends on uncertain parameters 
that involve the superconducting transition temperature 
and the entrainment factor. Here, all these uncertainties 
will be put into a single parameter e. In this case, e p can 
be written as: 



£ p — £p Pl 6.6el2- 

For typical neutron star numbers 
eSUx 10 6 ^, 



(4.4) 



(4.5) 



though its exact value is uncertain. Values for the vis- 
cosities are also needed. The electron and ordinary-fluid 
viscosities are given by |35| 



1.35 x 10 



19 



2.73 x 10 



14 



2 rji — 2 

Pl.5el4 J 8 > 



9/4 T -2 
F1.5el4 J 10 ■ 



(4.6) 



(4.7) 



Note that we define pi.5ei4 = p/(l-5 x 10 g ■ cm -3 ), 
T 8 = T/(10 8 K) (and so on for other numeric subscripts 
thoughout the rest of this paper). 

Finally, one important simplification must be pointed 
out. As discussed at the end of the last section, it is 
convenient to choose B — B r — constant. As already 
explained, this is unrealistic but should give the same 
qualitative results as a realistic field (see Mendell Q). 

Before finding the results for the critical angular ve- 
locity, it is easy to show that magnetic effects will be 
important for the superfluid model for the typical fields 
found in LMXBs. If the magn etic terms are larger than 
the viscous terms in Eq. (3.11) then magnetic effects on 
the MVBL length-scales will be important. Thus, the 
condition for this is 



V£v > Ve 



Pp 



(4.8) 



Using Eq. ([3/i]) in Eq. (L8) yields the following lower 



bound on the magnetic field, such that magnetic effects 
dominate the MVBL properties: 

B > (3.8 x 10 9 G)Kn 60 ^pl 5eU T- 2 p-l 6el2 ell e6 , (4.9) 

Thus, in the superfluid model radial magnetic fields ~ 
10 9 G (10 s K/T) 2 , not untypical of the fields of LMXBs, 
can affect the MVBL. Though, as we will show, the effects 
of mutual friction tend to counteract the magnetic effects. 
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FIG. 1. Temperature dependence of the critical angular ve- 
locities for ordinary-fluid neutron stars for B — 0, 10 9 (indis- 
tinguishable from B = 0), 10 10 , 10 11 , and 10 12 Gauss, from 
bottom to top. 



For the physical parameters and restrictions just dis- 
cussed, we solve Eq. for the critical angular velocity 
il c . The temperature dependence of the critical angu- 
lar velocity is presented in Figs. In these figures, 
the horizontal dashed line corresponds to 2f2 Q /3, which 
is the approximate maximum angular velocity for which 
mass-shedding occurs. Thus, the regions above this line 
are unphysical, and are kept only to illustrate the de- 
pendence of the curves for the range of magnetic field 
magnitudes typically found in neutron stars. Note that 
we extend the ordinary fluid curves into the superfluid re- 
gion (T < 10 9 K) and vice versa, given the uncertainties 
in the superfluid transition temperature. 

Figure l| shows the temperature dependence of the crit- 
ical angular velocities for ordinary-fluid neutron stars for 
B = 0, 10 9 (indistinguishable from B = 0), 10 10 , 10 u , 
and 10 12 Gauss, from bottom to top. Overall, it is seen 
that magnetic fields increase the dissipation rate, raising 
the critical angular velocity, and flatten out the curves. 
The slip factor, S, was set equal to 1 for all these curves. 
Smaller values of S lower these curves while preserving 
their shape. 

Figures || and || show results for the superfluid model 
for the respective cases when just the neutron vortices 
are pinned or when just the proton vortices are pinned. 
The dotted lines are for B = 0, 10 9 , 10 10 , 10 11 , and 10 12 
Gauss, from bottom to top when mutual friction is ig- 
nored (B n — 0), while the solid curves that diverge from 
these include the effects of mutual friction. The slip fac- 
tor, S, was set equal to 1 for all these curves. Smaller val- 
ues of S lower these curves while preserving their shape. 
Thus, while many of these curves are far above the break- 
up velocity of the star, they can be moved into the phys- 
ically allowed region by choosing a small value for S, 
which corresponds to weak pinning of the vortices. As 
in the ordinary-fluid model, magnetic fields increase the 




FIG. 2. Temperature dependence of the critical angular ve- 
locities for superfluid neutron stars with pinned neutron vor- 
tices for 5 = 0, 10 9 , 10 10 , 10 11 , and 10 12 Gauss, for dotted 
lines from bottom to top. Note that the dotted lines corre- 
spond to no mutual friction, while solid curves that diverge 
from these include mutual friction. 




FIG. 3. Temperature dependence of the critical angular ve- 
locities for superfluid neutron stars with pinned proton vor- 
tices for 5 = 0, 10 9 , 10 10 , 10 11 , and 10 12 Gauss, for dotted 
lines from bottom to top. Note that the dotted lines corre- 
spond to no mutual friction, while solid curves that diverge 
from these include mutual friction. 
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FIG. 4. Integrand that gives X for T = 10 9 K, B = 10 10 G, 
including mutual friction (lower curve) and excluding mutual 
friction (upper curve). 



dissipation rate, raising the critical angular velocity, and 
flatten out the curves. However, it is also seen that mu- 
tual friction tends to counteract the magnetic effects for 
high temperatures or high fields. 

The results shown here w ere c omputing numerically 
using t v as given by Eq. ( 3.32 ). However, we get 
the same resul ts with in a few tenths of a percent us- 
ing Eqs. ( 3.37 )-( j3~39 ). For our models, the constants 
in these equations are qoF — 4.37 x 10 _16 s _1 • G~ 2 , 
s" 1 / 2 , D OF = 1-56 x lO^s" 1 / 2 , 
G-\ while C SF = 0.261s- 1 / 2 , 
Asf — 0.0150s 1 / 2 , for pinned neu- 



Cof = 1.41 x 10- 6 *- 1 / 2 



;OF 
QSF 
D SF 



0.243 s 



'1/2 



tron vortices, and Csf = 4.31 x 10~ 5 s -1 / 2 , Dsf = 
1.05 x 10 _4 s -1 / 2 , for pinned proton vortices (we did not 
fit our data to dete rmine Asf fo r this case). Using these 
numbers and Eqs. (3.37)-(3.39), analytic expressions for 
t v for large B are 



33.9 s S~ Z B 



2 D -1 
12 ' 



(4.10) 



for the ordinarly-fluid model, and 
t v = 0.141 sS~ 2 B~ 1/2 



+4.85 x io- 5 sS- 2 TlBg /2 {n /n), 



(4.11) 



for the supefluid model when neutron vor tices are pinned. 
Substituting these equations into Eq. (4.1) and using 
Eq. (4.2) gives the following analytic expressions for the 
critical angular velocities 



(4.12) 



for the ordinarly-fluid model, while for the superfluid 
model with pinned neutron vortices we get 



^ = I.IQS^B 1 ' 12 . 

fin 



(4.13) 



when the first term on the right in Eq.(4.11) dominates, 
and 



^ = 9.7AS 2 /%- 2/5 B 9 1/w . 



(4.14) 



when the second term on the right in Eq.(4.11) domi- 
nates. 

We now explain why our results scale with B, T, and f2 
as they do. Note that in both the ordinary-fluid and su- 
perfluid models that the critical angular velocity becomes 
temperature independent for high temperatures (exclud- 
ing mutual friction). However, one surprising result is 
that mutual friction actually lowers the damping rate. 
To understand the various cases we need to study the 
characteristic length-scales associated with the boundary 
layer. 

The length-scales in the MBVL in the ordinary-fluid 
case have already been examined in Mendell f2~]}| . Here 
we will consider the superfluid case. The characteristic 
length-scales of the MVBL are determined by examining 
Eq. ( 3.11 ) and th e defi nitions of the length-scales given in 
Eqs. (3.35) and (3.36). Dropping subscripts for the pur- 
poses of discussion, since the MVBL damping rate is con- 
trolled by terms of the form i]\k\ 2 d = ij(4:Tt/X 2 + l/d 2 )d, 
it would seem that the smallest length-scales would de- 
termine the damping rate. However, as found by Mendell 
|pl|| , the ratio d/X turns out to be more important when 
magnetic effects are important. Also, only real frequen- 
cies, kQ, will be considered. This is valid either when the 
system is driven at a real frequency, or when the imagi- 
nary part of the frequency is small and can be ignored. 
(This is true for our models except near certain special 
angles, as will be explained.) Also, to further simplify 
the discussion, it will alwa ys be assumed that k > 0. 

First, consider Eq. (|3.1lD when B is less than the lower 
bounds given in Eq. ( [4.91) . The viscous term dominates, 
and the boundary layer thicknesses and wavelengths are 
given to lowest order by 



A 



1 



A =21 = _ / 2r l e 

71 2tt \K ± \Vn Pp 



(4.15) 



For typical neutron star numbers these length-scales are 



JUL 

2n 



47cm T ,_i n -i/2 -i/2 , R , 

• K | Pl.5el4-t 8 "6007rPp,6.6el2- (4.1b) 



The effect of mutual friction is to only add small cor- 
rections to the above, and thus does not matter in this 
case. This equation agrees in form with the standard 
result, but differs from previous results by numerical fac- 
tors of order unity due to the scaling with | K ± \ , and more 
importantly, because the proper scaling with the proton 
mass density is given here for the superfluid case, for the 
first time. This increases the VBL thickness, and con- 
sequently the damping time, compared to previous esti- 
mates for the superfluid case, e.g., by a factor of about 4.5 
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compared to the calculation given by Lindblom, Owen, 
and Ushomirsky [f20| . 

Mutual friction does matter for the case when B is 



greater than the lower bounds given in Eq. (4.E). How- 
ever, we first consider the case when mutual friction can 
be ignored. In this case K± is either purely real or purely 
imagi nary depending on the angle 9. Taylor expanding 
Eqs. (3.11) for this case the wave number is given by 



k± = K± i 



Un 2 



2 p p V^y 



(4.17) 



For the case of real K± the boundary layer length-scales 
to lowest order are 



1 



^cv 



A i 



2tt |^ ± | k !/2 ' 



(4.18) 



which is basically the distance a cyclotron-vortex wave 
travels in one rotation, and 



VcvPp 



ds/ri \k±\k?/ 2 n 2 Ve 



(4.19) 



The subscripts indicate whether these quanties depend 
on purely magnetic, or a ratio of magnetic to viscous 
quantities. Substituting in values for the parameters 
gives 



As 
2tt 



33 cm i /2 x 

o nl L, 



,1/2 



\K±\k>-/*' 



3.8e9"6007r c 1.4e6' 



(4.20) 



65 cm 3/2 2 2 3/2 Pp,6.6el2 , . , 

a Bh ~ \ tF. l„3/2 - D 3.8e9 s '6007r J 8 £ 1.4e6~ ' ^- Z1 J 



|Jf±|K 3 / a ' 



Pl.5el4 



For the case of imaginary K± the roles of A and d become 
interchanged. Thus, for this case, the results are as in 
Eqs. (pQ^)-(|42l|) with 



d 
A 



A_ 

2nd. 



(4.22) 



It is apparent that for large B that the MVBL damping 
rate is dominated by terms of the type 



B 1 ' 2 . 



(4.23) 



This explains the regions where the damping rate be- 
comes independent of the temperature and angular ve- 
locity, and the scaling of t v with _B -1 / 2 in the first term 
on the right side Eq. ( 1.11 ). Physically, the reason that 
magnetic fields increase the damping rate is that while 
magnetic forces increase the size of the boundary layer 
many wavelengths of cyclotron-vortex waves fit into the 
boundary layer thickness (A < d) and there is a large 
amount of shear associated with these waves. (Similar 



explanations are given in Mendell pll] for the ordinary- 
fluid case.) 

Attention will now be given to cases tha t occ ur near 
certain special angles in reference to Eq. ( 3.12 ). The 
angles are defined by 



k ± 2cosA± = 0, 
k ± 27Cosy± = 0, 



Pnl + Pp 



k ± 27CosZ± = 0. 



(4.24) 
(4.25) 

(4.26) 



The latter two of these angles are important to under- 
standing the effects of mutual friction (and mutual fric- 
tion smooths out singular behavior near these angles). 
For the models used in this paper, these angles are lo- 
cated at: X- = 1.23, X+ = 1.91, Y- = 1.39, Y+ = 1.74, 
Z_ = 1.24, and Z + = 1.90. 

First, for near X±, the small imaginary part of k 
becomes important (in both the ordinary-fluid and su- 
perfluid models). As explained in Mendell [ pl[ , to un- 
derstand this, make the replacement k — > k + i/(flr), 
and note that fir ^> 1 fo r the situations of interest. It is 
seen that k± oc y/i/ (Or) and that the angular derivative 
dk±/dQ becomes large near X± (infinite at X± when the 
imaginary part of k is ignored). However, as shown in 
Mendell these regions in the ordinary-fluid model 
make little contribution to the VBL damping rate. In 
the superfluid model, the proximity of X± to Z± tends 
to narrow the singular behavior near X± . Thus, correct- 
ing for the small imaginary part of k near 9 = X± is 
not done in this paper, nor has it been made in previous 
studies. 

Now consider the case when mutual friction is in- 
cluded. Its effects are easiest to understand near the 
angles Y± and Z±, which are special only in the super- 
fluid model. The imaginary part of k could also be im- 
portant near these angles, but note that mutual friction 
effects are more important as long as r > l/(2jflB n )- 
Thus, mutual friction effects dominates near Y± and Z± 
when the damping time is greater than the lower bound 
t > 2.7s/(7n6007r<6 nj i.0e-4)' It is found that this lower 
bound is valid for the results presented in this paper. 
Also, recall that mutual friction is only important for the 
case B is greater than the lower bound given in Eq. ([4.9]). 
With these caveats in mind, first consider the effects of 
mutual friction when 9 is near Y±. One can show that 
the characteristic MVBL length-scales are 



A_B/g„ _ j2 / p n pp Vgy 



which for typical neutron star numbers is 
d 



(4.27) 



l B/B n 



Ag/ B „ _ 670 cm x/2 _i 1/2 

~ 2-K ~ 7k 3 - 8e9 6007r £ 1.4e6 



1/2 1/2 
Pn,1.434el4Pp,6.6el2 w-1/2 
X ~ t 'n,1.0e-4- 

Pl.5el4 



(4.28) 
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Next, for 6 
order are 



Z±, the MVBL length-scales to lowest 



A 



BB„ 

2ir 



(7-1) 




2p 2 B n V cv 



PnPp 



(4.29) 



which for typical neutron star numbers is 



IB-S„ 



2tt 



(3.2cm)7 

(T-l)^ 2 !!» 1Ae6 



Pl.5el4 



■ 1/2 1/2 
Pn,1.434el4/ 3 p,6.6el2 



B 



1/2 



(4.30) 



We see that the length-scales again depend on the dis- 
tance a cyclotron- vortex wave travels in one rotation, but 
that d — \/2ir. Thus we do not have many wavelengths 
within the boundary layer, but only one wavelenth. We 
see that the length-scales near Z± are similar to the case 
when mutual friction is ignored. However, near Y± the 
length-scales are much larger than in the case when mu- 
tual friction is ignored. This reduces the dissipation rate 
and is the main reason that mutual friction reduces the 
critical angular velocities as shown in Figures || and ||. 
Also note, that when d — X/2tt the dissipation varies as 
r]\k\ 2 d ~ r)/d, and so the dam ping time varies as d/rj. 
For d in Eqs. (4.28) and ( 4. 30] ) we see that this implies 
t v ~ i? 1 / 2 T 2 /fl. Thus, we have explained the scaling 
of the mutual friction corr ectio n term with B, T, a nd 
foun d em pirically in Eq. ( |3.39| ) and used in Eqs. ( 4.11 ) 
and ( 4.14 ). The overall effect of mutual friction is shown 
in Fig. |J, which again shows that mutual friction lowers 
the dissipation rate, and smooths out singular behavior 
near the special angles considered here. 



V. CYCLIC EVOLUTION OF ACCRETING 
NEUTRON STARS 

We now turn our attention to the differential equa- 
tions governing the spin cycles of accreting neutron stars. 
In the context of the r-modes, these cycles have been 
discussed by Levin [03[ , Anderson^t al. [JiTj, Wagoner, 
Hennawi, and Liu ||34|| , and Heyl |35|. Previously Wag- 
oner |3l|] suggested that the spins of these stars could 
be limited by gravitational-radiation from an unstable 
mode, while Bildsten proposed an alternative model 
in which the gravitational-radiation is generated due to 
an accretion-induced asymmetry of the star. 

First, as in Owen et al. we define the dimensionless 
r-mode amplitude a as the maximum perturbed fluid 
velocity at the equator of the star divided by HR. (Note 
that we do make corrections in this section for the fact 
the r-mode is confined to just the core of the star but 
that this does not change the definition of a.) Second, 
following the analysis of Wagoner, Hennawi, and Liu ]34| ] 
(which builds on the work of Owen et al. Iq] and Levin 
p3[), the following equations govern the evolution of the 



r-mode amplitude, a, the angular velcocity, f2, and the 
temperature T: 



1 da 
a dt 



— F — F 

— r g 1 v 



F g K c a 2 -^F a , 



1 dfl 
Q~dt 



-2F g K c a 2 + F a , 



dT 

C(T)— = W diss 



K n Mc 2 — L„(T). 



(5.1) 



(5.2) 



(5.3) 



These equations, of course, are highly simplified in that 
they ignore differential rotation, assume infinite heat con- 
duction, and so forth, but give a qualitative overview of 
how the evolution should proceed. 

We now discuss the various terms in these equations. 
Heating due to dissipation is given by 



Wdiss = JMR 2 n 2 a 2 F v . 



(5.4) 



The constant J is defined in Eq. (3.4) of Owen et al. || 
and relates to the canonical angular momentum of the 
r-modes. Here we recompute this constant for r-modcs 
confined to just the core of the star and find its value is 
J = 0.01255. The constant K c , called Q in Owen et al. 
Q and defined below their Eq. (3.4), is related to J and 
the moment of inertia the star (which is unchanged even 
when the r-modes are confined to the core). However, 
since J changes, K c also changes, and we find K c = 0.072. 
Next, note that Wagoner, Hennawi, and Liu |34|] define 
1 — Kj as the fraction of the canonical r-mode angular 
momentum that contributes to the total physical angu- 
lar momentum of the star. In these equations we have 
adopted the case that the canonical angular momentum 
of the r-modes contributes zero physical angular momen- 
tum to the star, i.e., Kj = 1. We find in our case, as Wag- 
oner, Hennawi, and Liu [ pMf find for the cases they study, 
that the exact value of Kj is not important to the results 
since the r-modes saturate quickly and no more physi- 
cal angular momentum can go into them, if any at all 
does to begin with. Furthermore, Levin and Ushomirsky 
p6| show that the canonical angular momentum of the 
r-modes contributes zero physical angular momentum to 
the star for a toy model, further justifying our use of 
Kj = 1. Next, note that the above equations are valid 
only for small a. This is OK even when the mode satu- 
rates as long as the saturation amplitude is small. Several 
studies suggest the saturation amplitude is indeed small 
p6| , p7| . However, for comparison with earlier work we do 
consider a few cases where the saturation amplitude is 
a = 1 , and for these few case only we instead use Eq. (5) 
in Levin |33| to evolve £1 when a = 1. As in Wagoner, 
Hennawi, and Liu |54|, but using the MBVL damping 
rate which dominates for our models, we set 
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FIG. 5. Spin frequency vs. temperature evolution for ordi- 
nary-fluid neutron stars for B = G and B = 10 10 G (both 
for S = 1.0), and for B = 10 11 G (for S = 0.1). The curves 
are explained in the text. 



where tqr is given by Eq. (42), while we set 
F =1 
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FIG. 6. Spin frequency vs. temperature evolution for su- 
perfluid neutron stars with pinned neutron vortices and for 
B = G and S = 0.01. The curves are explained in the text. 



(5.6) 



where t v is the MVBL damping time given by Eq. (3.32). 
The inverse accretion time-scale, F a , is related to the 
accretion rate, M, by 



F a = 



1 



M 



5 x 10 6 yr J \ lO- 8 M /yr 



(5.7) 



where the observed accretion rates in LMXBs are in the 
range 10- n M Q /yr < M < 10" 8 M Q /yr |57|. We take the 
accretion rate to be constant. Even if the accretion rate 
were not constant, while the r-modes are saturated the 
accretion rate has little effect on the evolution. Note that 
heating due to accretion is proportional to the constant 
K n = lx 1(T 3 Q. Finally, in Eq. the heat capacity 
and neutrino cooling rates are given by H||,[34],^7j , 



C= 1.47 x 10 38 erg • K _1 T 8 , 



•8 



T ' T 16 
±J brcm- L S > 



for the ordinary-fluid case, and 

C = 7 x 10 36 erg • KT 1 ^ 

^brcm 1 8 ' 



T — J' T 6 



(5.8) 
(5.9) 



(5.10) 
(5.11) 



for the superfiuid case. The constants L' URCA and L' brem 
are given by |57J i 



L'urCA = 1-5 x 10 32 ergs \ 
L' hiem = 1.4 x lO^ergs- 1 . 



(5.12) 
(5.13) 



The results of evolving Eqs. (|5.l|)-(|5.3|) are the cyclic 
curves shown in Figs. |^-^[ The monitonically decreas- 
ing smooth curve through the middle of each cycle is 
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FIG. 7. Spin frequency vs. temperature evolution for su- 
perfiuid neutron stars with pinned neutron vortices and for 
B — 10 9 G and S — 0.005. The curves are explained in the 
text. 
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FIG. 8. Spin frequency vs. temperature evolution for su- 
perfiuid neutron stars with pinned neutron vortices and for 
B = 10 1 G and S = 0.01. The curves are explained in the 
text. 

the critical angular velocity curve. As time increases, a 
star progresses clockwise around the cycle. Note that 
for each cycle there is a maximum and minimum angu- 
lar velocity and a maximum and minimum temperature. 
At the maximum angular velocity the r-mode becomes 
unstable at the minimum temperature, i.e., we define 
Sl max = fl c (T m - ln ). At the minimum angular velocity the 
r-mode becomes stable at the maximum temperature, 
i.e., we define fl mnl = r2 c (T max ). We define the change 
in angular velocity by f2 max — ^min = Af2. We will de- 
scribe some cycles as "thin" if Af2 <C f2 m in, otherwise 
we will describe the cycle as "fat" if fi m i n -C f^max- For 
each cycle the star goes through roughly 4 distinct stages. 
Each stage corresponds to one side of the roughly quadru- 
lateral curve that corresponds to one cycle ]5S| ]. These 
stages have been previously discussed by Levin p3| and 
Anderson, et al. [^7| (see also Heyl p5|). Here we extend 
their results to include the superfluid case, magnetic ef- 
fects, and thin cycles. First we give typical numbers for 
our models. After this we discuss the figures in detail. 

The first stage corresponds to the vertical left side of 
the cycle, which corresponds to the star spinning up by 
accretion with heating due to accretion balanced by cool- 
ing due to neutrino emission K n Mc 2 = L„(T). Thus, the 
minimum temperture is related to the accretion rate by 



T mill = 2.8 x 10 s K 



M 



1/8 



lO- 8 M /yr 



for the ordinary-fluid case, and 



T 



.6 x 10 8 K 



M 



1/6 



10- 8 M Q /yr 



(5.14) 



(5.15) 



for the superfluid case. The time-scale for accretion to 
spin up the star is given by solving Eq.( |5.2| ) when F a 
dominates: 



+ - 1 1 n » 

tspmup — m 

r a i l Ti 



5 x 10 6 yr 



lO- 8 M /yr 
M 



In 



fir 



(5.16) 



It will be useful to note that for thin cycles the spinup 
time is approximately, 



5 x 10 6 yr 



M 



An 

"min 



(5.17) 



The second stage corresponds to the the top of the 
cycle. The beginning of this stage occurs when the r- 
mode becomes unstable and the mode amplitude grows 
to a saturation amplitude a sa t- This happens on the 
gravitational growth time-scale of 



1 , Qsat 

— m 



9.3 hr 



600vr/s 

ax 



In 



Osat 



(5.18) 



As is done in Wagoner, Hennawi, and Liu |34|], we assume 
a residual r-mode amplitude is continuously excited by 
stochastic processes and do not allow a to evolve be- 
low a m in = 10~ 12 . As the mode grows, viscous dissipa- 
tion heats the star to a maximum temperture for which 
JMR 2 Q 2 a1 at F v = L V (T). Thus, the maximum temper- 
ature is determined by the saturation amplitude. If we 
overestimate F v by replacing it with ^(fimax) we find 



T max < 5.0 x 10"' 



V10- 4 / 



V6007t/e 



for the ordinary-fluid case, and 



T ma ,x < 1.9 x 10 9 K 



VlO- 4 / li 



4/3 



(5.19) 



(5.20) 



.10-V V600tt/s, 
for the superfluid case. Next, we can solve Eq. (5.3) when 



the term with F v dominates to get the characteristic time 
for dissipation to heat the star 



the 



1 



2W dil 



[C(T a 



C(T min )T min )] , (5.21) 



If we make the approximations T max ^ T m ; n and F v < 
-Fg(^max), estimates for this time-scale are 



or ( as 
9oyr — 

J V 10 



3/2 / O 



.10- 4 / V600tt/s 
for the ordinary-fluid case, and 



the 



69 yr 



/_Osat_\ 

V10- 4 / 



-4/3 / O 



V600tt/s 



-16/3 



(5.22) 



(5.23) 



for the superfluid case. Note that this time-scale domi- 
nates the top of the evolution cycle. 
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The third stage occurs when angular-momentum radi- 
ated away a s gr avitational radiation spins down the star. 
Solving Eq. ( ^2 ) when the term with F g dominates gives 



spindown 



1~GR 



a. 



For thin cycles, 



spindown 



7.4 x 10 5 yr 



An /io 



2 / 600tt/ s n 6 

V " m i t 



(5.24) 



(5.25) 



while for fat cycles 

^spindown 



This stage dominates the right side of the evolution cy- 
cles. Note that stages 2 and 3 begin at the same time. 
Because they have very different durations they can be 
studied separately. 

The fourth stage begins when the star spins down to 
the critical angular velocity curve. The amplitude of the 
mode then decays away on the dissipation timescale 



amp 



1 ln QWt 

Fv ^min 

„„, /600tt/s 
9.3 hr ' 

\ * 'min 



In 



a sa t 



(5.27) 



where we have used F v ~ ^(fimin) during the damping 
to obtain the numerical estimate for this time-scale. The 
last stage is dominated by the time it takes the star to 
cool back to T m i n . Solving Eq. (5.3) when the L v domi- 
nates gives 



(10 s K) 
6Z( 



URCA 



C(T n 



C(T n 



7^7 

max 



4.9 x 10 5 yr 



10 8 K 

-*min 



(5.28) 



for the ordinary-fluid case, and 

(10 8 K) 6 TC(T min ) C(T n 



cool 



4L1 
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7^5 

max 



4.0 x 10 fa yr 



10 8 K 

-L mm 



(5.29) 



for the superfluid case. This stage dominates the bot- 
tom of the evolution curves. Note that stage 4 really 
coincides with stage 1 of the next cycle. Because they 
have somewhat different durations they can be studied 
separately. 

Thus, as Levin 33 1, Anderson, et al. ]17[] , and Heyl 
p5[ conclude, we see that the time that the r-mode 
spends saturated (and thus actively radiating potentially 



detectable gravitational waves) is dominated by the spin- 
down time-scale, while the time that the r-mode spends 
with minimum amplitude (and thus radiating insignifi- 
cantly) is dominated by the spinup time. Thus, defining 
r as the fraction of LMXBs in the active phase gives 



r = t. 



spindo wn /^s pinup ; 



as in Levin 



using Eqs. ( gig ) and djTgg ), 



For fat cycles, 



< .072 



10 



x 2 /500tt/s s 6 



■ 



M 

lO" 8 M /yr 



, (5.30) 



where In(f2 max /f2 m i n ) was treate d as a fact or of order 
unity. For thin cycles, using Eqs. (5.17) and (5.25), 



< .44 



io- 



Osat 



M 



/500tt/s n 

V fimin J llO- 8 M Q /yr 



(5.31) 



The first of these ratios is the same as what Levin |33| ] 
found, though scaled differently here. Note that for these 
equations we have scaled fi min with the observed mini- 
mum value 250 Hz, rather than a typical value of 300 Hz 
to get upper bounds on r. We see that LMXBs spend 
about 6 times as much time radiating if their evolution 
is controlled by thin rather than fat cycles. However, 
current observations suggest that the mostly thinly clus- 
tered set of LMXBs probably have f2 m i n /27r < 250 Hz and 
/2tt > 350 Hzj|§30|M. 



This case falls between the 



a 

cases that Eqs. ( 5.31 ) a nd (5.30 ) ar e vali d. Instead, using 
these numbers in Eqs. (5.16) and (5.24) gives r = 0.18. 
(It is possible of course that a subset of these might be- 
long to an even more thinly clustered set of LMXBs; on 
the other hand, some pulsars have Q,/2n > 600 Hz in- 
dicating that some LMXBs have been able to spin up 
well beyond 350 Hz.) Levin |53| gives an estimate of 10- 
100 as the number of strongly accreting LMXBs in our 
galaxy. This would indicate that if the r-modes instabil- 
ity controls the cycles of LMXBs there is a good chance 
some LMXBs in the galaxy are currently radiating as 
long as the saturation amplitude is small. Levin |35| l and 
Anderson, et al. ]l7[ | come to more pessimistic conclu- 
sions, based on a larger saturation amplitude. However, 
we should not be too optimistic either, and treat our 
values for r as upper bounds. Many LMXBs have accre- 
tion rates far lower than 10 _8 Mo/yr, lowering the real 
chances that any LMXB currently radiates. Decreasing 
the size of the saturation amplitude increases the odds, 
but also makes the r-modes harder to detect. We discuss 
detectability again at the very end of this section. (Heyl 
p5| , at the same time this paper was written, indepen- 
dently considered small saturation amplitudes. That pa- 
per estimates the number of LMXBs in the galaxy based 
on the observed number of millisecond pulsars and comes 
to a similar but perhaps more optimistic conclusion than 
presented here.) 

We now discuss the cycles in each figure. In Fig. ^ 
curves are shown for the ordinary-fluid case. The dashed 
curves with large cycles are for M = 10 _8 Af Q yr _1 , 
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"sat = 1, and 5 = 1. The largest of these is for B = G 
and the slightly smaller cycle is for B = 10 10 G. The 
lower dashed curve showing a thin cycle is for 5 = 0.1 
and B = 10 11 G. This is the largest buried field that could 
exist in the crust of an LMXB [B7| , and we choose the slip 
factor to put this cycle within the observed range of spin 
frequencies of LMXBs. The thin solid lines are cycles for 
the same parameters as the dashed cycles, except for a 
more realistic saturation amplitude of a sa t = 10 -4 . Re- 
cent studies indicate the saturation amplitude is small, 



However, 
the entire 



and perhaps much smaller than 10~ 4 ]26],|27 
when the saturation amplitude is this smal 
cycle moves into the superfluid region. Since the tran- 
sition temperature is uncertain we do not exclude these 
curves as a possibility. The main conclusion drawn here 
is that only a very large fields can produce thin cycles in 
the ordinary-fluid case. 

In Figs. HH results are shown respectively for the su- 
perfluid model with pinned neutron vortices for B = 0, 
10 9 , and 10 10 Gauss, and for S = 0.01, S = 0.005, and 
S = 0.01. In each case, the slip factor is choose to put 
the curves within the observed range of spin frequen- 
cies observed in LMXBs. The case of pinned proton 
vortices would result in similar plots, but for different 
values of the slip factor. The dashed curves correspond 
to M = lO- 8 M yr _1 and a sat = 0.1. The thick solid 
curves are for the same accretion rate, but a sat = 10~ 4 . 
The thin solid curves are for M = 10~ 11 M Q yr~ 1 and 
«sat = 10~ 4 . The surprising main conclusion drawn here 
is that B ~ 10 9 has the largest effect on narrowing the 
cycles. This is because this value of B tends to flatten 
the critical anglular velocity curve in the relevant tem- 
perature range, while mutual friction causes it to vary 
with temperature again for larger B fields. 

In general, lowering the saturation amplitude moves 
the lower right-hand corner of cycle to the left along the 
critical angular velocity curve, while increasing the ac- 
cretion rate moves the upper left-hand corner of a cycle 
to the right along a critical angular velocity curve. Thus 
adjusting these can make a cycle more or less narrow, as 
can the magnitude of the magnetic field. This makes it 
hard to use current observations of LMXBs to place con- 
traints on the internal physics of these stars. However, 
if a class of LMXBs with similar internal magnetic fields 
could be identified, one prediction is that weakly accret- 
ing LMXBs should exhibit a wider range of spins than 
strongly accreting LMXBs. In any case, magnetic fields 
can make the cycles thinner, and are thus potentially im- 
portant to the explanation of the observed clustering of 
LMXB spin frequencies. Thus, if the appropriate mag- 
netic fields are present in a class of LMXBs it will tend 
to make the cycles thinner and increase the number of 
currently radiating sources. 

Given the uncertainties we cannot predict with any 
certainty if any LMXBs are currently radiating. How- 
ever, no matter the odds, the detection of gravita- 
tional radiation from the r-modes is not ruled out 
for gravitational-wave detectors currently coming online 



(with their planned enhancements) as long as the satu- 
ration amplitude is not too small. Using Eq. (4.9) from 
Owen et al. Q but correcting for the fact that in our 
models the r-modes are confined to the core, the average 
dimensionless gravitational-wave amplitude for a source 
a distance D away is, 



h = 7.6 x 10" 



/ «sat \ / Q V 
VlO" 4 / V6007r/ S y 



/lOKpc 



(5.32) 



For this fudicial value of h, and assuming an integrated 
observation time of 10 7 s, the optimal signal-to-noise 
(SNR) ratios are in the range 0.57 - 4.4 for the LIGO 
noise curves given in Owen et al. |^j. (Since it is not 
possible to keep the detector optimally aligned with the 
source actual SNRs would be 5 — 10 times smaller. How- 
ever, advanced detectors can greatly improve the SNR in 
narrow frequency bands.) 

Thus, assume that gravitational radiation is detected 
from an LMXB in our galaxy. The radiation can be iden- 
tified as due to the r-mode instability if the observed fre- 
quency is in the correct ratio to the spin frequency of 
the LMXB. The measured gravitational-wave amplitude 
then determines the r-mode saturation amplitude (if the 
distance to the LMXB is known) which determines the 
right side of a spin cycle. If the accetion rate is known 
the left side of a spin cycle is known. It is unlikely to ob- 
serve an LMXB along the top or bottom of a spin cycle. 
However, the maximum (or minimum) spin frequencies, 
which give the top (or bottom) of a spin cycle must be 
larger( or smaller) than the observed values for these. 
Thus, it might be possible to place limits on the interior 
magnetic fields of LMXBs with known accretion rates de- 
pending on how well the allowed range of spin frequencies 
can be determined. 



VI. CONCLUSIONS 

We have studied the r-modes in accreting neutron stars 
with magneto-viscous boundary layers for ordinary-fluid 
and superfluid models. The ordinary-fluid model is dis- 
cussed in detail in Mcndell |2l]]. We have shown, within 
our approximations, that no solution to the MHD equa- 
tions exists when both the neutron and proton vortices 
are pinned. However, solutions can exist for cases when 
just one species of vortex is pinned. In these cases, the 
instability of the r-modes to gravitational-wave emission 
can limit the spins of neutron stars (though strongly 
pinned neutron vortices would completely suppress the 
instability). If both the neutron and proton vortices 
are completely unpinned, the results would be basic ally 
that found in Lindblom and Mendell |Q and Levin [^3| . 
In general, magnetic fields increase the dissipation rate 
and flatten the critical angular velocity vs. tempera- 
ture curves in both models. However, mutual friction in 
the superfluid model tends to counteract the magnetic 
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effects for high temperatures or high fields. Further- 
more, if the r-mode instability controls the spin cycles 
of LMXBs we have shown several things. As in Levin 
p3| Anderson, et al. |T^| ; and Heyl [p5f , we show that 
decreasing the saturation amplitude greatly increases the 
odds that at least one LMXB in our galaxy is currently 
radiating gravitational waves. We also show that mak- 
ing the spin cycles thinner further increases the odds by 
up to a factor of 6 compared to previous estimates. Fi- 
nally, we show for the first time that magnetic fields, by 
flattening the critical angular velocity curves, make the 
spin cycles of LMXBs thinner and that this increases the 
fraction of time an LMXB spends radiating gravitational 
waves for B > 10 11 G in the ordinarly-fluid case and for 
B ~ 10 9 G in the superfluid case. Previous studies pred- 
icated that probably no LMXBs are currently radiating 
in our galaxy, based on large values of the saturation 
amplitude a sa t- However, the subsequent studies of Wu, 
Matzner, and Arras ||(| and Arras et al. predict that 
a will be small. In any case, the previous studies agree 
with our result that the fraction of currently radiating 
LMXBs approaches one if a sat < 10~ 4 and if at least 
10 strongly accreting LMXBs exist in our galaxy. How- 
ever, if a sa t is much larger than 10 -4 then probably no 
LMXBs are currently radiating; if a sa t is much smaller 
than 1CP 4 the r-modes may be undetectable. (Detection 
of a source in our galaxy is difficult, but not impossible for 
a as small as 10~ 4 .) However, if gravitational waves are 
detected from the r-modes of an accreting neutron star, 
it might be possible to place limits on the interior mag- 
netic fields of LMXBs with known accretion rates based 
on their observed allowed range of spin frequencies. 

Of course many caveats must be added to our conclu- 
sions. Our evolution graphs show that the fluids may 
go through a transition from superfluid to ordinary-fluid 
during the heating phase, and vice versa during the cool- 
ing phase. Vortices would probably not have time to 
migrate out of the core during the spindown phase, and 
thus superfluid effects could remain important above the 
transition temperature; but the details of what happens 
to the boundary layer during such a phase transition 
could become very messy. This is only one complica- 
tion that we have ignored. The neutrons could be in 
the superfluid phase while the protons are in the nor- 
mal phase (for example see Easson and Pethick [ p9| ) and 
vice versa. Another possibility is that the protons could 
form a type I superconductor in the intermediate state 
rather than a type II superconductor in the vortex state 
H . We have also only studied the effects of mutual fric- 
tion for one value of the entrainment factor; special val- 
ues exist that would completely change our results fl48| . 
Furthermore, we have ignored interactions between neu- 
tron and proton vortices (see j40|) and other dissipative 
effects, such as hyperon bulk viscosity [ p2|j2^j2^ , |25|| . Fi- 
nally, nonlinear effects are ignored, such as the winding 
of magnetic fields lines, which could be important espe- 
cially for B > 10 10 G Qn]]6|]. The winding of field 
lines would introduce another temperature independent 



damping mechanism, and thus probably would comple- 
ment our results. Clearly, an exact understanding of the 
r-mode instability is very complicated. 

Finally, we remark on one final interesting question: 
If the r-modes control the spins of LMXBs, could the 
observation of gravitational waves from some of these 
objects distinguish between ordinary-fluid neutron stars, 
superfluid neutron stars, and strange stars? The recent 
study of Andersson, Jones, and Kokkotas [[33| suggests 
that strange stars in LMXBs may emit persistent gravi- 
tational waves from saturated r-modes. We have shown 
that neutron stars in LMXBs, in agreement with Levin 
p3[ , Anderson, et al. jl?]], and Heyl pq ], would emit grav- 
itational waves for only part of their spin cycle. If we do 
not observe gravitational radiation from any LMXBs we 
probably do not learn much, since so many mechanisms 
can suppress the instability. However, if the r-modes are 
observed from some LMXBs, and the current theoretical 
understanding holds up, then the fraction of LMXBs that 
radiate would distinguish between neutron and strange 
stars. A detailed comparison between theory and obser- 
vation would be needed to distinguish between ordinary- 
fluid and superfluid neutron stars, however. Obviously, 
further work is needed to understand the various possi- 
bilities. 
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APPENDIX: SUPERFLUID MVBL SOLUTION 

We will now solve the corrective differential equations 
governing MVBL dissipation in a superfluid neutron star 
with arbitrary B-field and without our ignoring a factor 
of cos 2 in Eq. (3.4). For such a model, we obtain the 
system of differential equations: 



A5v e + BS^ = FSi + G5v 
-BSv e + ASv^ = JS'v + K5v 4 

C5v° + DSv* = FSv + GSv 
-DSv e + ESv* = JS'v + KSv 



HSw + ISw 
LSw + MSw 4 
How + low 
LSw + MSw 4 



where 



A = ipn 

B = -2p cos 9 

C = ippK + 2'fpB n cos 2 8 
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ippK + 2"/pB 

1 [ -ip p V c 



2 

cv 



B< 
~~B 



{B B f 



B 2 



G = 


iPpVgy 


K n 2 


H = 


P 


I = 


Pn n 
— 7G 

P 


J = 


G 


K = 


1 f-i, 



{B r ) 2 B e B<> 



B 4 



ip p V c 



2 

cv 



L = 



M 



Pn T 
— 7J 

P 



B< 
~B 



2 r 



1 - 



{B*f 



B 2 



As in the simplified case, this system of differential equa- 
tions only has two linearly independent variables: 



Sv 8 



aSw 9 + b5w* 
c8w e + ddw* 



where 



b = 



d = 



AC 


+ BD 


A 2 


+ B 2 


AD 


- BE 


A 2 


+ B 2 


BC 


-AD 


A 2 


+ B 2 


AE 


+ BD 


A 2 


+ B 2 



This leaves us with two linearly independent systems 
of differential equations. Assuming the same exponential 
behavior for 5v and 5w as in the simplified model, we 
find 



f3 ± yfp 2 - 4a 7 



2a 



where 



= PS- QR 
= PE + SC 
= D 2 + CE 



DQ - DR 



and 



P = aF + cG + H 
G = bF + dG + I 
R = a J + cK + L 
S = bJ + dK + M. 
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